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In a three-body system of identical bosons interacting through a large S-wave scattering length a, 
there are several sets of Efimov features related by discrete scale invariance. Effective field theory was 
recently used to derive universal relations between these Efimov features that include the first-order 
correction due to a nonzero effective range rs. We reveal a simple pattern in these range corrections 
that had not been previously identified. The pattern is explained by the renormalization group for 
the effective field theory, which implies that the Efimov three-body parameter runs logarithmically 
with the momentum scale at a rate proportional to rs/a. The running Efimov parameter also 
explains the empirical observation that range corrections can be largely taken into account by 
shifting the Efimov parameter by an adjustable parameter divided by a. The accuracy of universal 
relations that include first-order range corrections is verified by comparing with various theoretical 
calculations using models with nonzero range. 
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Introduction. Dramatic experimental and theoretical 
progress in few-body physics has been stimulated by the 
realization that Feshbach resonances can be used to con¬ 
trol the strength of interatomic interactions in ultracold 
atomic gases. The low-energy universality that arises 
when scattering lengths are much larger than the range 
of interactions is of particularly broad interest. It implies 
that for systems as disparate as atoms, hadrons, and nu¬ 
clei, dimensionless combinations of few-body observables 
are the same, despite orders of magnitude differences in 
the length and energy scales. 

A particularly fascinating class of low-energy univer¬ 
sal behavior is Efimov physies, which is characterized by 
discrete scale invarianee [1, 2]. The ability to control 
the strength of interatomic interactions has allowed the 
observations of various aspects of Efimov physics in ultra¬ 
cold atoms [3] . The simplest example of Efimov physics is 
the Efimov effect: in the unitary limit where the scatter¬ 
ing length a is infinite, there are infinitely many three- 
body bound states with an accumulation point at the 
3-particle scattering threshold [4] . These Efimov trimers 
have binding energies ET,n whose limiting behavior is 

ET,n —^ nl/m as n ^ 00 , (1) 

where m is the particle mass and A is called the discrete 
sealing faetor. The limit in Eq. (1) defines a three-body 
interaction parameter /t*. In the case of identical bosons, 
A = = 22.6944 with sq = 1.00624. The discrete 

spectrum in the unitary limit in Eq. (1) reveals that low- 
energy physics in the three-particle sector depends not 
only on the scattering length a but also on /c*. Efimov 


physics can also be revealed through discrete scale in¬ 
variance at finite values of a. For example, the negative 
scattering lengths a-^n at which the Efimov trimers cross 
the three-boson threshold have the limiting behavior 

O-^n — >0-X"k~^ as n —>• 00 , (2) 

where = —1.50763 [5] is a universal number. 

The discrete scale invariance that characterizes Efimov 
physics becomes exact, with the limits in Eqs. (1) and (2) 
replaced by equalities, only in the limit of zero-range in¬ 
teractions. In that limit, few-body observables depend 
only on a and k*. The interactions in real physical sys¬ 
tems always have nonzero range. Experiments naturally 
involve the deepest Efimov trimers, for which range cor¬ 
rections are largest. In order to make quantitative ex¬ 
perimental tests of Efimov physics, it is important to 
understand the range corrections in detail. 

Range corrections can be studied theoretically by cal¬ 
culating Efimov features in various models for inter¬ 
atomic interactions and comparing them with the univer¬ 
sal zero-range predictions. Kievsky and Gattobigio dis¬ 
covered empirically that range corrections can be largely 
taken into account by making substitutions for a and k* 
in zero-range formulas [6] . Their substitution for the Efi¬ 
mov parameter is k* —^ -I- T/a, where T is determined 
empirically for each observable and each system. 

A systematic theoretical approach to the problem is to 
organize range corrections into an expansion in powers 
of the range tq [7]. In the two-body sector, the only new 
parameter that enters through second order in ro is the 
S-wave effective range rs- Bedaque et al. showed that 
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three-body observables do not depend on any additional 
three-body parameter at first order in rg [8], but they 
made the implicit assumption that the scattering length 
a is fixed. Ji, Phillips, and Platter (JPP) showed that if 
variations in the scattering length are considered, there 
is an additional three-body parameter at first order in 
To [9]. They developed an effective field theory (EFT) 
framework for calculating range corrections [9, 10], and 
used it to derive universal relations between Efimov fea¬ 
tures, such as ET,n and a_,n, that are accurate to next- 
to-leading order (NLO) in rg/a. 

In this paper, we reveal a pattern in the first-order 
range corrections to Efimov features that was not rec¬ 
ognized in Refs. [9, 10]. We demonstrate that the pat¬ 
tern has a simple renormalization-group interpretation 
in terms of a running Efimov parameter that runs with 
the momentum scale at a rate proportional to rs/a. The 
empirical shift of the Efimov parameter used in Ref. [6] 
to incorporate range corrections into zero-range formulas 
can be identified as the expansion of the running Efimov 
parameter to first order in Vg/a. We compare the pre¬ 
dictions of the NLO universal relations with results for 
Efimov features in specific models. 

Efimov Features. The Efimov trimers can be labelled 
by an integer n. In the 3-atom sector, the most dra¬ 
matic Efimov features associated with the n**' branch of 
Efimov trimers are (a) the binding momentum kt,™ = 
of the Efimov trimer in the unitary limit 
a = ±oo, (b) the negative scattering length a_^„ at which 
the Efimov trimer crosses the 3-atom threshold, (c) the 
positive scattering length at which the Efimov trimer 
disappears through the atom-dimer threshold, and (d) 
the positive scattering length at which there is an 
interference minimum in the three-body recombination 
rate at threshold. 

The three-body interaction parameter k* defined by 
Eq. (1) is approximately equal to the binding momen¬ 
tum of the Efimov trimer labelled n = 0. The binding 
momenta are Kr.n = in the zero-range limit. JPP 

showed that this equation remains exact, with no range 
corrections, at first order in rg/a [11]. In the universal 
zero-range limit, the ratio of any pair of Efimov features 
is a universal number. The leading order (LO) universal 
relations are 


^i,n — A , 
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expressed using Eq. (3) as l/a_^g = OZ^K^+Iulrs, which 
defines a nonuniversal number I. JPP showed that the 
range expansion to first order in the range for any other 
Efimov feature can then be expressed as 

l/az,n = A -t- {fi^n T : (4) 

where and ? 7 i_„ are universal numbers. JPP calcu¬ 
lated many such numbers to at least three digits [9, 10]. 

There is a pattern to the dependence of the universal 
numbers in Eq. (4) on the number n labelling the branch 
of Efimov trimers that was not identified in Ref. [10]. 
The range expansion in Eq. (4) can be expressed in the 
much simpler form 

ai^n = + {Ji - ni7)rs, (5) 


where a = 1.095 is a universal number. The differences 
between the coefficients Ji are also universal numbers: 
J+ —J_ = 0.548, J* —J_ = 1.250. We will refer to Eq. (5) 
as the NLO range expansion for the Efimov feature. 

Renormalization. The n in the range correction in 
Eq. (5) can be understood as a logarithmic dependence 
on the scale at which the observable Oi^n is measured, 
suggesting a renormalization-group interpretation. We 
therefore discuss the renormalization of the EFT used 
to calculate range corrections in Refs. [9, 10]. The La- 
grangian density for that EFT, which has an atom field 
■0 and a molecule field is 

—^ {iffipij) -\- h.c) -|- h(j)^if. (6) 

v2 

The parameters g and A can be tuned as functions of the 
ultraviolet cutoff, A, so that a and Vg have the desired 
values. In the three-body sector, the ultraviolet cutoff 
can be implemented as an upper limit A on the loop mo¬ 
mentum in a modified Skorniakov-Ter-Martirosian equa¬ 
tion [12] for the scattering of three bosons [13]. The 
three-body coupling constant can be expressed as h = 
2mg^E[/A^, where H is a dimensionless log-periodic func¬ 
tion of A. Its dependence on A in the zero-range limit 
was determined up to a numerical factor in Ref. [13]: 


= 0.879 ""[*° ‘°] , 

sm[sg log(A/A,) — arctansg] 


( 7 ) 


where 6»_ = -1.50763, 9+ = \9_\/VX = 0.316473 [5], and 

= 0.0707645 [1]. 

First-order Range Corrections. JPP developed an 
EFT framework for calculating range corrections as strict 
expansions in powers of rg [9, 10]. At first order in rg, 
two Efimov features are required as inputs, with at most 
one being a trimer binding momentum (or k*). A simple 
choice for the two Efimov features is k* and a__g. The 
deviation from the zero-range prediction for a__g can be 


where A* = 0.548 k*. The multiplicative constant in 
Eq. (7) was first determined in Ref. [14]. The dependence 
of on A at first order in the range has the form [10] 


HiA) 


Ho{A/kC) + /iig(A/K*)Ars 
7i7o(A//c0 log(A//rg) -h hii{A/nf) 



where po is a momentum scale and Hq = {Ad/dA) Flo is 
the logarithmic derivative of i7g in Eq. (7). We will not 
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need the analytic forms of the log-periodic functions /iio 
and hii. JPP gave an analytic expression for 7 which is 
within 5% of the value we find numerically, 7 = 0.351. In 
Ref. [10], ^0 was chosen so that a second Efimov feature 
in addition to k*, such as a_.o, had no first-order range 
corrections. The order of magnitude of the required value 
of ^0 is then l/|a__o|- 

Running Efimov Parameter. The term proportional to 
log(A/^o) in Eq. ( 8 ) represents a logarithmic violation 
of discrete scale invariance. It can be absorbed into the 
zero-range coupling constant by changing the argument 
of iJo to The power of ^0 in the 

argument can be canceled by replacing k* by a running 
Efimov parameter defined by 

K,(/io,a) = K* . (9) 

Unnecessarily large logarithms in the range corrections 
to an observable dominated by the momentum scale Q 
can be avoided by expressing the observable in terms of 
R^,{Q,a), a, and instead of k*, a, and r^. 

In the NLO range expansion in Eq. (5), the term lin¬ 
ear in n gives unnecessarily large range corrections if n 
is large. Natural scaling variables in the zero-range limit 
are the inverse scattering length 1 /a and the energy vari¬ 
able K = sig-D.{E){m\E\/. The momentum scale 
Q for an observable with energy E at scattering length 
a is Q = {K^ + 1/a^)^/^. The momentum scales for 
the Efimov features a_,„, a+^„, and are l/|a__„|, 
l/a+_„, and •\/2/a*,n, respectively. Since a factor of 
is not essential, these can be summarized via the LO 
universal relation as Q ~ If we replace k* in 

the LO universal relation in Eq. (3) by k;*(( 5, a) where 
a = X^OiKf^ and expand in powers of rg, we obtain 
a term —n'yrg log A. This matches the —narg term in 
Eq. (5) provided a = 7 log A. Our numerically deter¬ 
mined values of 7 and a satisfy this condition to within 
numerical accuracy. Replacing k* in the LO universal re¬ 
lation by A"0iKL^) therefore includes the 

NLO correction proportional to nr^. 

Similarly, renormalization-group improvement of the 
NLO range expansion can be obtained by eliminating At* 
in Eq. (5) in favor of the appropriate running Efimov 
parameter and demanding agreement to first order in r ^: 

At*-! ^ (JO) 

where Ji = Ji + ylogj^ij. The differences between the 
coefficients Ji are universal numbers: J+ — J_ = 0.000, 
— J- = 0.177. We will refer to Eq. (10) as the RG- 
improved NLO range expansion for the Efimov feature. 
For n = 1, the NLO range expansion in Eq. (5) and 
the RG-improved NLO range expansion in Eq. (10) have 
the same parametric accuracy: higher-order range cor¬ 
rections are suppressed by a factor of (At*rg)^. For large 
n, the NLO range expansion has corrections of order 


n^(A-"At*rs)^. In the RG-improved NLO range expan¬ 
sion, those higher-order range corrections that are en¬ 
hanced by a factor of n for every factor of r^ are summed 
up to all orders. Thus the higher-order range corrections 
to Eq. (10) are at most of order n{X~^K^rsY■ 

Comparisons with Models. Successive Efimov features 
have been calculated in several models with nonzero 
range. We can use the results to illustrate the accu¬ 
racy of the NLO range expansion in Eq. (5) and the RG- 
improved NLO range expansion in Eq. (10). The NLO 
range expansion can be used to express ai_„+i/(Aa„_i) as 
a linear combination of K^rg and JiK^^rg with universal 
coefficients. For any three ratios of successive Efimov 
features, there is a linear combination with universal co¬ 
efficients in which A«*rs is eliminated. We refer to this 
equation as an NLO universal relation. If two of the ra¬ 
tios are taken as inputs, the third can be predicted to 
NLO accuracy without knowing rg. Similarly the RG- 
improved NLO range expansion can be used to derive 
a universal relation that expresses log(ai_„_|_i/(Aai^ji)) as 
a linear combination of two other such expressions with 
universal coefficients. 

TABLE I: Ratios of successive Efimov features a_,n+i/a-,n 
divided by the discrete scaling factor A. The ratios of the 
features calculated by Deltuva [15] are compared to the pre¬ 
dictions of Eq. (5) (NLO) and Eq. (10) (RG-NLO) using the 
numbers in square brackets as inputs. 



( (1 — ^ 7J,-(-1 / d — ^ 71) / A 

n 

0 

12 3 

4 

Ref. [15] 

0.7822 

0.9665 0.9976 0.9999 

1.0000 

NLO 

[0.7822] 

[0.9665] 0.9975 0.9998 

1.0000 

RG-NLO 

[0.7822] 

[0.9665] 0.9975 0.9998 

1.0000 


In Ref. [15], Deltuva calculated the scattering lengths 
at which universal tetramers cross the 4-boson thresh¬ 
old for identical bosons interacting through a separable 
Gaussian potential. He also gave accurate results for the 
ratios of the scattering lengths at which 6 

successive Efimov trimers cross the three-boson thresh¬ 
old. These ratios divided by A are given in Table 1. They 
rapidly approach 1 as n increases. In Table I, we have 
taken the ratios a__„+i/a__„ for n = 0 and 1 as inputs, 
and then used NLO universal relations to predict the ra¬ 
tios for n = 2, 3, and 4. The predictions are in excellent 
agreement with the results calculated by Deltuva. The 
RG-improved NLO universal relation gives the same pre¬ 
dictions to four digits. 

Schmidt et al. calculated multiple Efimov features for 
a three-parameter model in which the interaction is a 
transition between an atom pair and a molecule with a 
Gaussian form factor [16]. The three independent pa¬ 
rameters can be specified by a, r^, and a dimensionless 
parameter Sres = obtained by dividing the range 
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TABLE II: Ratios of successive Efimov features ai^n+i!ai^r. 
divided by the discrete scaling factor A. The ratios calcu¬ 
lated by Schmidt et al. [16] are compared to the predictions 
of Eq. (5) (NLO) and Eq. (10) (RG-NLO) using the numbers 
in square brackets as inputs. The three blocks are for models 
with Srea = 100, 1, aud 0.1. 



(«-, 

.1+1 ja-^ 

n)/A 


)A 

n 

0 

1 

2 

0 

1 

2 

Ref. [16] 

0.753 

0.962 

0.998 

0.175 

1.764 

1.029 

NLO 

[0.755] 

[0.962] 

0.997 

-8.814 

1.150 

1.032 

RG-NLO 

[0.755] 

[0.962] 

0.997 

0.0002 

1.206 

1.034 

Ref. [16] 

1.008 

0.998 

0.9998 

0.757 

0.983 

1.001 

NLO 

[1.008] 

[0.998] 

0.9998 

-0.431 

0.986 

1.002 

RG-NLO 

[1.008] 

[0.998] 

0.9998 

0.240 

0.986 

1.002 

Ref. [16] 

1.156 

1.012 

1.0007 

1.188 

0.938 

0.991 

NLO 

[1.156] 

[1.012] 

1.0007 

0.449 

0.869 

0.990 

RG-NLO 

[1.156] 

[1.012] 

1.0008 

0.916 

0.887 

0.990 


To of the Gaussian form factor by a length r* determined 
by its strength. They presented results for three sets of 
parameters with Sres = 100, 1, and 0.1. Their results 
for a-^n+i/0'-,n and divided by the discrete 

scaling factor A are given in Table II. We have taken the 
ratios a-s/a-^ and a_^ 2 /a-,i as inputs, and then used 
NLO universal relations to predict the other ratios. All 
predictions for n = 2 are in good agreement with the 
results of Ref. [16]. The predictions for n = 1 are also 
in reasonable agreement, except for the case Sres = 100. 
The large error in this case can be attributed to the fea¬ 
ture a_.o being outside the window of universality. The 
RG-improved NLO universal relation gives slightly better 
predictions for n = 1 and dramatically better predictions 
for a* 4 /a*_ 0 j for which the NLO universal relation yields 
unphysical negative values in the Sres = 100 and 1 mod¬ 
els. 

Shift in Three-body Parameter. In Ref. [6], Kievsky 
and Gattobigio calculated the binding energies of Efi¬ 
mov trimers and the atom-dimer scattering length in two 
models with a Gaussian two-body potential and with or 
without a Gaussian three-body potential. They made 
the empirical observation that range corrections could 
be largely taken into account by making substitutions in 
the zero-range formulas, which are functions of a and k* 
only: (a) replace a by the inverse binding momentum of 
the universal dimer or virtual state, (b) replace the Efi¬ 
mov parameter k* by k* -I- T/a, where T is a parameter 
that is determined empirically for each observable and 
each system. In Ref. [17], Garrido et al. showed that this 
prescription also works for the three-body recombination 
rate at threshold. The accuracy of the prescription was 
verified only in limited regions of a and for a few specific 
models. The miming Efimov parameter provides a the¬ 


oretical justification for the prescription. The expansion 
of the running Efimov parameter in Eq. (9) to first order 
in the range is 

K,(Q,a) Ri K, [1 - 7log(Q/K*)rs/a]. (11) 

If the slow logarithmic dependence of the momentum 
scale Q on a is ignored, this has the same form -\-T/a 
as the empirical shift in the Efimov parameter introduced 
by Kievsky and Gattobigio. 

The NLO range expansion can be used to identify uni¬ 
versal relations between the empirical constants T for dif¬ 
ferent 3-body observables. The prescription in Ref. [17] 
for Efimov features associated with the second Efimov 
trimer is o^qAtT,! TTi = 9i. According to the NLO range 
expansion in Eq. (5), the ratio (T^ — Tj)/{KT.irs) should 
be the universal number Jj — Ji. The predicted universal 
value of J, — J+ is 0.702, while the two models considered 
in Ref. [17] give the results 0.671 and 0.753. We consider 
this to be reasonable agreement. The corresponding ra¬ 
tios involving r_ display large discrepancies between the 
two models, and do not agree as well with the universal 
predictions. 

Experiment. By measuring the first Efimov feature 
at different Feshbach resonances in ^^^Gs atoms and com¬ 
paring with previous measurements with other atoms, 
the Innsbruck group discovered a correlation to within 
20 % between and the coefficient of the van der 
Waals tail —Ce/r® of the interatomic potential; a__o ~ 
—OcvdW: where rvdw = [18]. This corre¬ 

lation, which can be called van der Waals universality, 
was subsequently verified theoretically [19, 20]. Van der 
Waals nniversality gives a sharp prediction for the effec¬ 
tive range: Ts = 2.79rvdw [21]. It thus predicts a narrow 
range of values for the coefficient J_ in Eq. (5). 

To use the NLO range expansion to predict Efimov fea¬ 
tures for a specific bosonic atom, the required inputs are 
the scattering fength a, the effective range rg, and two 
measured Efimov features. There are not many atoms for 
which there are accurate measurements of three or more 
Efimov features. Severaf Efimov features have been ob¬ 
served for both ^Li and ^^^Cs, but there are complications 
in ^Li from signihcant variation of and in ^^^Gs from 
multiple Feshbach resonances. Comparison of NLO pre¬ 
dictions with experimentally measured Efimov features 
will be presented elsewhere [22]. 

Summary. We have shown that universal relations 
which account for a nonzero effective range are able to 
predict Ehmov features with higher accuracy. A simple 
pattern in the NLO range expansion can be interpreted 
in terms of a running Efimov parameter. The running 
Efimov parameter also explains the empirical findings of 
Gattobigio et al. that range corrections can be incorpo¬ 
rated into universal zero-range results through a shift in 
the Efimov three-body parameter. 
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